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The exact solvability of several nuclear models with non-degenerate single-particle energies is 
outlined and leads to a generalization of integrable Richardson-Gaudin models, like the sw(2)-based 
fermion pairing, to any simple Lie algebra. As an example, the so(5) ~ sp(4) model of T = 1 pairing 
is discussed and illustrated for the case of Ge with non-degenerate single-particle energies. 
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Exactly solvable models (ESM) built upon a dynami- 
cal symmetry have a long history of providing important 
insights into the structure of nuclei. The two main advan- 
tages of ESM are: (1) They can describe in an analytical 
or exact numerical way a wide variety of elementary phe- 
nomena. (2) They can be and have been used as a testing 
ground for various many-body approaches. The simplest 
example of ESM is the rank-one su{2) model of fermions 
in one orbit or in several degenerate orbits with a con- 
stant pairing interaction, often used to introduce nuclear 
superconductivity (see e.g., Ref. Q). 

In general, a quantum system has a dynamical symme- 
try if the Hamiltonian can be expressed in terms of the 
Casimir operators of a chain of nested algebras. A typical 
example of a model with rank-two Lie algebra dynamical 
symmetry is Elliott's sm(3) model, introduced to describe 
the phenomenon of nuclear deformation j^]. Elliott's 
Hamiltonian is a linear combination of the quadratic 
Casimir operator of the su(3) Lie algebra (involving a 
quadrupole-quadrupole interaction) and of the Casimir 
of its so(3) subalgebra, associated with angular momen- 
tum. Lie algebras with higher rank lead to more complex 
ESM like the so(5) model of T = 1, isovector pairing Q, 
the so(8) model of T = 0, 1 isoscalar and isovector pair- 
ing HJali Ginocchio's model with so(8) and sp{6) struc- 
ture [g , also known as the Fermion Dynamical Symmetry 
Model (FDSM) Q. The three dynamical symmetries of 
the Interacting Boson Model (IBM) Q provide another 
example of this approach. 

The concept of quantum integrability goes beyond the 
limits of the dynamical symmetry approach. A quantum 
system is integrable if there exist as many commuting 
Hermitian operators (integrals of motion) as quantum 
degrees of freedom -Sij . The set of Casimir operators of a 
chain of nested algebras fulfills this condition. There are, 
however, also well-known examples of integrable models 
without dynamical symmetry |lClj |. 

Usually dynamical symmetry models are defined for 
degenerate single-particle levels. Changes in the single- 



particle energies break the dynamical symmetry but may 
still preserve integrability. The pairing model with non- 
degenerate single-particle levels, with an exact solution 
found by Richardson in the sixties, represents a unique 
example of an ESM with such characteristics jll|. Re- 
cently, Richardson's model has been shown to be inte- 
grable by finding the complete set of integrals of mo- 
tion constructed in terms of the generators of the su(2) 
algebra [Tjl- Subsequently, more general exactly solv- 
able pairing models, both for fermions and for bosons, 
called the Richardson-Gaudin (RG) models, have been 
proposed [l3- Since then, a great deal of work has been 
devoted to understand the properties of these ESM and 
to apply them to a wide variety of problems in nuclear, 
condensed- matter, atomic, and molecular physics |l4| . 

The aim of this Letter is to present the generalization 
of the RG models to those symmetry algebras that have 
given rise to the well-known ESM in nuclear physics. In 
most cases this generalization will allow non-degenerate 
single-particle energies, as well as other symmetry break- 
ing one-body operators. As an example, we shall dis- 
cuss the exact solution of the so(5) proton-neutron (pn) 
isovector pairing model with non-degenerate orbits, for 
which an exact solution has been proposed by Richard- 
son [13 • However, it has been recently shown that the 
Richardson's solution was incorrect [Ig ■ 

We also mention that so(5) has been proposed as 
the symmetry underlying high Tc-superconductivity |l7| . 
The exactly solvable RG models discussed in this Let- 
ter could be used to generalize so(5) condensed-matter 
models jig by the explicit addition of non-degenerate 
single-particle symmetry-breaking terms. 

We begin by introducing a set of commuting operators, 
the RG operators (integrals of motion) Ri 19]: 
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where the index i [j) refers to the i-th (j'-th) copy of a 
Lie algebra C with generators Xf , Zi are fixed parameters 
which later will be related to the single-particle energies, 
and ^ is a generic element of the Cartan subalgebra of £ 
(see below). The scalar product is defined through the 
>C-invariant metric tensor gap — c^pCp^ where c^p are 



the structure constants of £, e.g., Xi ■ Xj = X^g^fsX' . 
Since the R^ commute, any function of these operators 
can be used as a model Hamiltonian for an integrable 
system. In particular, any linear combination of the Ri 
is integrable and is at most quadratic in the generators. 
For any simple algebra C, the eigenvalues of the Ri 
have been given in Ref. 1 191 o r can be derived from the 
Gaudin-algebra approach |2(]| |: 
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where the parameters Cnn are solutions of the generalized 
Richardson equations [I^ |23| 
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The primed sum means that the singular term (a, a) = 
(6, P) is omitted and L is the number of copies of the 
algebra £ of rank r. The Cartan subalgebra £^ C £ 
has the elements hi, . . . ,h^ and the p" are the compo- 
nents of ^ in this generally non-orthonormal basis. The 
7r° are the components of the r simple roots ttI, . . . , tt'' 
of/: in the Cartan- Weyl basis, [hs,E''] = tt'^^E''. Each 
Ai is a vector of highest weight (usually of the fun- 
damental representation) for the i-th copy of C and 
the expression A^ • tt" corresponds to the eigenvalue of 
iJ" = J2s ""s ^'' in this highest-weight state of the i-th 
copy of £. The p° are, in a sense, the strength com- 
ponents of the symmetry-breaking one-body operator ^. 
Finally, the M" are positive numbers related to the eigen- 
values m" of H'^ in the eigenstate of the integrals of mo- 
tion d), M" = J2^iA^ ■ T^" - <)■ 

Equations Q and © were derived by Asorey et al. 
under the assumption that C is the Lie algebra of a sim- 
ple, simply-connected, compact Lie group. We arrived at 
the same equations within Ushvcridizc's framework |20j 
which assumes that £ is a singular semi-simple Lie al- 
gebra (i.e., a classical algebra or one of the exceptional 
algebras Eg or E-j). Since Eqs. jSJl and |j3Jl involve a 
dot product, well defined for any Lie algebra, one may 
expect that the equations themselves are valid for any 
semi-simple Lie algebra and possibly for an even larger 
class. In fact, it is desirable to find a formulation of these 
equations that involves arbitrarily chosen Cartan algebra 
generators instead of the simple roots of the algebra. 

Although one can use any function of the Ri as a 
Hamiltonian, the following particular linear combination 
yields a simple expression for the eigenvalues: 
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With the same linear combination of the eigenvalues 
ri |j2Jl and using the Richardson equations lj3Jl, one can 
show that the eigenvalues of Q are linear in the spectral 



parameters Ca^a with coefficients 5°' that depend on the 
symmetry-breaking strengths p". 

The Hamiltonian Q) shows how the global C- 
symmetry represented by the total second-order Casimir 
operator C2 is broken by the first term containing the 
elements S,i of the Cartan subalgebras £,p. In the general 
case of Zi ^ Zj the first term does not commute with C2 . 
Thus, the eigenstates of (0)) are spread over different irre- 
ducible representations of C and they are not eigenstates 
of the total Casimir operator C2 . 

Using Cartan's classification of semi-simple Lie alge- 
bras, one can now generalize many nuclear physics mod- 
els within this framework (see TableQJ . Note that models 



TABLE I: Algebras associated with some nuclear physics 
models; ~ denotes isomorphisms. R&R is the symplectic 
model of Rosensteel and Rowe l2l|. *FDSM Lie algebras ^. 



rank n 


su{n + 1) 


so(2n + 1) 


sp{2n) 


so(2n) 


1 


pairing 


~ su{2) 


~ su{2) 


~u(l) 


2 


Elliott 


T — 1 pairing 


~ so(5) 


~ su{2) e su{2) 


3 


Wigner 


soiiy 


R&R* 


~ su{A)* 


4 


su(5) 


so(9) 


sp(8) 


T — Q,l pairing* 



based on a fermion realization of the type X^p ~ 'Ax'^a^ 
such as sp{2n) or so{n), are well adapted for the gener- 
alization towards non-degenerate single-particle energies 
because in that case the Cartan generators are sums of 
fermion number operators. In contrast, models based on 
realizations of the form a^^ap, such as su{n), contain dif- 
ferences of number operators. 

To illustrate the present formalism, we shall discuss 
the example of T = 1 pairing in systems with protons 
and neutrons and non-degenerate single-particle levels. 
The dynamical symmetry limit of pn-pairing in a single 
degenerate shell was first considered by Flowers 3] as an 
important extension of the original pairing problem be- 
tween identical particles. It yields a seniority classifica- 
tion of protons and neutrons in jj coupling in terms of pp, 
nn, and pn pairs. In doing so, the concepts of seniority v 
(the number of nucleons not in pairs coupled to angular 
momentum zero) and of reduced isospin t (the isospin of 
these nucleons) are established. The application of this 
so(5) formalism has given rise to many analytic reduction 
formulae for shell-model matrix elements |23| . 

For simplicity, in the RG extension of the so(5) dy- 
namical symmetry model we will restrict ourselves to se- 
niority V = (and consequently reduced isospin t — 0), 
though states with broken nucleon pairs can be easily in- 
corporated in the formalism. We begin by specifying the 
Cartan generators of so(5): 
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where p labels protons (p) or neutrons (n), i is the i-th 
copy of so(5) and i is the time reversal of the state i. We 
identify hi with the third component of the isospin and 
/i2 with the nucleon number. 

The positive-root vectors arc the T ~ \ creation opera- 
tors {h^_^^^ = njn|,6j_, = {n\p\+p\n\)lV2,b\^^^ ^ p\p\} 
plus the isospin-raising operator T+_,; = {p\ni+p\ni) / ^/2. 
The simple-root vectors are {h_n,T+^i} and h_i ^ is the 
singular-root vector 12y|. These, together with the conju- 
gate operators, close the so(5) algebra. The isospin sub- 
algebra sut(2) C so(5) is generated by the so(5) elements 
that are not related to b_^ ^. That is, sut{2) is generated 
by T+^i, T-^i = {T+.j)^ , and [T+,i,T_,i] = Tq^j = /ii,.;. 

Since so(5) is a rank-two algebra, there are two types of 
spectral parameters: ei^a and 62^/3 which shall be denoted 
as Wa and ep. The upper bounds Ml and M2 for the 
indices a and /?, respectively, are related to the isospin 
T and the total number of pairs N via the expressions 
Ml = N -T and M2 = N. The scalar products of the 
simple roots are 7r2-7r2 = 2, ttI -ttI = 1, and 7r2-7rl = —1. 
In the spherical shell model protons and neutrons occupy 
single-particle states with quantum numbers (j, m). The 
index i then corresponds to jm and i to jm. Alterna- 
tively, due to the rotational symmetry, the angular mo- 
mentum j can be used as a label instead of i but then 
the corresponding degeneracy Vtj = (2j -I- l)/2 should be 
taken into account. Finally, the weights A; are the same 
for all i and are those of the fundamental representation 
of so(5), that is, A-7rl = and A-7r2 — 1. Inserting these 
in Eq. Q, together with the choice pi = 0, p2 = ^l/ff, 
and Zi = 2£i, we obtain the generalized Richardson equa- 
tions for the T = 1 pn-pairing: 
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The particular case of flj — 1 was derived by Links et 
al. [23 using the algebraic Bethe ansatz. Each solution 
of the equations |SJ| gives an eigenstate of the pn-pairing 
Hamiltonian 



i/p„ = V2£,(/j2„ 
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with eigenvalues E = X)a=i ^a- The spectral parameters 
e are interpreted as pair energies as in the case of su{2) 
pairing. However, due to the larger rank of the so(5), a 
new set of spectral parameters wg appears in the equa- 
tions (O. These new parameters w are associated with 
the su{2)t isospin subalgebra. For each possible isospin 
T there are N — T new w-parameters. The meaning of 
this new set of parameters w, which do not appear in the 
expression for the eigenvalues, becomes evident when an- 
alyzing the eigenstates of ©. The Bethe ansatz for the 



so(5) eigenstates of the RG model is a factorized product 
wave function. It consists, for A^ < ^ fi^, of a neutron- 
pair product state related to the lowest-weight state of 
sut(2) and a isospin-raising product operator state that 
guaranties the required total isospin of the ansatz: 

N-T N 

\w,e-e) = W r+K)[]feli(e„)|o), 

7=1 a— 1 

where the spectral dependence of the generators is given 
by the expressions 



^+(^'7) = E 
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Note that application of T+ on b_i results in 6q, and ap- 
plication of T_|_ on 6q results in bj^^ which provides us with 
Ushveridize's linear combination of powers of 6j, |20|. 

Further insight into the structure of the equations Q 
can be gained by using the classical electrostatic anal- 
ogy Hi- They can be considered as the equihbrium 
condition for a 2D classical electrostatic problem involv- 
ing three sets of charged particles: (a) the orbitons with 
charges — $7^/2 and fixed positions 2£j\ (b) the pairons 
with unit positive charges and free positions ec,; and (c) 
the new w particles with unit positive charges and free 
positions wg. In addition, there is a uniform electric 
field in the vertical direction (the real axes) with strength 
~ 1/(7. Furthermore, the w particles do not interact with 
the orbitons and do not feel the electric field, which can 
bee seen from the second equation in lO). 

As an example, we performed numerical calculations 
for ^"^Ge (4 protons and 4 neutrons) using the single- 
particle energies (in MeV) 63/2 = 0.00, £5/2 = 0.77, 
ei/2 = 1.11, and eg/2 = 3.00 and two pn-pairing 
strengths, g = 0.1 (weak) and 0.37 (strong). We first 
checked, in this non-trivial case, that the solutions of (Q 
are indeed the exact solution of the Hamiltonian (|SJ by 
direct comparison with a standard nuclear shell-model 
calculation. 

Figure ^ shows the solutions for the lowest T = 0, 
1, and 2 states. The T = solution corresponds to 
the ground state, while the T — 1 and T = 2 solu- 
tions are excited states in ^''Ge. As in the su{2) pairing 
case, the different configurations can be classified in the 
weak-coupling limit. As can be seen from Fig.^ at weak 
coupling four pairs occupy the P3/2 (the four pairons are 
close to the P3/2) in the T=0 state. This configuration is 
not allowed for T = 1 and T = 2 states due to the Pauli 
principle. Correspondingly, three pairons are close to the 
lowest P3/2 orbiton and the fourth pairon approaches the 
orbiton 75/2. In all cases the w particles are intertwined 
with the pairons to minimize the Coulomb interaction. 
The number of w particles {N — T), together with the 
initial configuration at weak coupling, defines each eigen- 
state of the pn-pairing Hamiltonian. As g increases, the 
free charges expand under the influence of the electric 
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FIG. 1: Two-dimensional representation of the e pairon and 
w particle positions in ''^Ge for the lowest-energy states with 
isospin r = 0, T = 1, and T = 2. The left panel corresponds 
to weak coupling g = 0.1 and the right panel to strong cou- 
pling g — 0.37. The black squares represent the two lowest 
orbitons {ps/2, fs/i) with positions (2£j), the grey circles are 
the e pairon positions, and the open circles are w particle 
positions. 



field and their mutual interaction. The solutions are sub- 
ject to numerical instabilities due to singularities arising 
when a real pair energy e crosses a single-particle energy 
or when real e and w parameters cross. An example of 
the first class of crossings can be observed in Fig. ^ for 



T = 2 where the pairon associated with the f^/2 orbiton 
at weak coupling goes down with increasing g and crosses 
the P3/2 orbiton. The T — 1 case shows an exchange of 
positions in the real axis of a pairon and a w particle as 
an example of the second class of singularities. The first 
class of singularities were already present in the su{2) 
pairing case and they precluded the practical use of the 
exact solution for a long time. Recently, a new method to 
overcome this numerical problem was proposed [23 ■ We 
beheve that the same procedure can be used to treat the 
second class of singularities as well, making feasible the 
exact solution of the so(5) model for very large systems. 

In summary, we have presented the generalization of 
RG models to arbitrary semi-simple Lie algebras, which 
include most of the dynamical-symmetry models of nu- 
clear physics. The generalized RG models allow the intro- 
duction of one-body symmetry breaking terms like non- 
degenerate single-particle energies. As an example of this 
approach, we gave the exact solution of the so(5) pn- 
pairing model. We emphasize that the exact solution for 
large systems with so(5) symmetry could be of great im- 
portance in condensed-matter physics in addressing the 
phenomenon of high Tc-superconductivity |l7l [ia| . Fi- 
nally, the treatment of higher-rank algebras like sp(6) 
and so(8) opens the possibility of exact nuclear struc- 
ture calculations with more realistic quantum integrable 
models. 
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